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ABSTRACT 

Aims. Gravitational settling of 22 Ne in cooling white dwarfs can affect the outcome of thermonuclear supernovae. We investigate how 
the supernova energetics and nucleosynthesis are in turn influenced by this process. We use realistic chemical profiles derived from 
state-of-the-art white dwarf cooling sequences. The cooling sequences provide a link between the white dwarf chemical structure and 
the age of the supernova progenitor system. 

Methods. The cooling sequence of a 1 Mq white dwarf was computed until freezing using an up-to-date stellar evolutionary code. 
We computed explosions of both Chandrasekhar mass and sub-Chandrasekhar mass white dwarfs, assuming spherical symmetry and 
neglecting convective mixing during the pre-supernova carbon simmering phase to maximize the effects of chemical separation. 
Results. Neither gravitational settling of 22 Ne nor chemical differentiation of 12 C and l6 have an appreciable impact on the properties 
of Type la supernovae, unless there is a direct dependence of the flame properties (density of transition from deflagration to detona- 
tion) on the chemical composition. At a fixed transition density, the maximum variation in the supernova magnitude obtained from 
progenitors of different ages is ~ 0.06 magnitudes, and even assuming an unrealistically large diffusion coefficient of 22 Ne it would be 
less than ~ 0.09 mag. However, if the transition density depends on the chemical composition (all other things being equal) the oldest 
SNIa can be as much as 0.4 magnitudes brighter than the youngest ones (in our models the age difference is 7.4 Gyr). In addition, our 
results show that 22 Ne sedimentation cannot be invoked to account for the formation of a central core of stable neutron-rich Fe-group 
nuclei in the ejecta of sub-Chandrasekhar models, as required by observations of Type la supernovae. 
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O i 1. Introduction stars (lHamuv et al.ll996ll2000tlWan? et al.l2006l) . This observa- 

y—i ' tional finding implies that the age difference of the progenitors 

O ! Thermonuclear explosions of accreting white dwarfs are thought is one of the origins of SNIa diversities. The age effect may have 

^ — 1 ■ to be the origin of several of t he most violent phenomena known: its origin in either the intrinsic dependence of the observational 

Type l a supernovae (SNIa) dHovle & Fowled! 19601 iBranch et alj properties of individual supernovae on the age of their progen- 
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1995; Hillebr andt & Niemeven 12000 1 : iRopkd l2006h . .Ia super- itors (this could be called a pure age effect) or the occurrence 
nova e (Bildsten et al. 2007), and even some Type lb supernovae of different channels leading to SNIa, each one characterized by 
dPerets et al.ll2010h. The discovery of the accelerated ex pansion slightly different properties and time histories (e.g. white dwarf 



of the Universe (Ri ess et al.ll9 98: Perlmutt er et all 19991) has put mergers versus single degenerate progenitors), or even in obser- 



SNIa at the center of attention of cosmologists, as a complete vational biases. In this work, we focus on a possible pure age 

knowledge of any systematics affecting SNIa luminosity is nec- effect caused by the chemical differentiation that occurs during 

essary to achieve the precision required to m easure the equation white dwarf cooling prior to the supernova event. 

of state of dark matter (see, e.g. Tonry 2005). ~, , • . c ff . 

\ & i ui v Observational constraints on the nature of the age effect 

One of the most intriguing systematics concerning SNIa come from statistical studies of th e properties of large nu mbers 

variability is the correlation suggested by the observations be- of SNIa in different environments. Gallagher et all d2008l) found 

tween supernova brightness and progenitor age, the so-called a strong correlation suggesting that SNIa in galaxies whose pop- 

"age effect" (there is another type of age effect that correlate s ulation is older than 5 Gyr are ~ 1 mag fainter at maximum 

the SNIa rate with age, see for instance Nom oto et alj 12000). than those found in galaxies with younger populations. They 

The luminosity of SNIa is related to the morphological types of concluded that the time since progenitor formation primarily de- 

their host galaxies: brighter SNIa tend to occur in spiral galaxies termines the production of Ni, the main radioactive nuclide 

with younger stellar populations, while most of the fainter events powering the light curve. A similar co nclusion was obtained by 

occur preferentially in early-type galaxies with relatively older iHowell et al]d2009h . lNeill et al.l d2009) pointed out what appears 
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to be a threshold average population age of 3 Gyr above which 
a host is less likely to produce SNIa with 56 Ni masses greater 
than ~ 0.5 Mq. Below this age threshold, there appears to be 
little correlation between 55 Ni mass and host age. These find- 
ings tightl y constrain the pos sible origin of the age effect. For 
instance, Krue ger et alj f2010) attributed the luminosity-age cor- 
relation to an increase in the central density of the white dwarf 
prior to the SNIa explosion, during the accretion phase of the 
progenitor system. However, they only explored cooling times 
be low ~ 1 Gyr, whic h are much less than the threshold suggested 
bv lNeill etalJd200l. 

iTimmes et all (l2003h and lTonrvl d2005l) argued that sedimen- 
tation of 22 Ne over a timescale of several (up to ~ 7 - 8) Gyr 
provides a natural time-dependent mechanism modulating the 
luminosity of a white dwarf explosion. The sequence of events 
can be divided into three phases. First, the primary component 
of a binary system evolves into a white dwarf. For the range of 
white dwarf progenitor masses, the timescale of the first phase is 
on the order of < 1 Gyr. Second, there is a variable period of time 
during which the secondary remains in the main-sequence and 
the white dwarf cools and acquires a (perhaps partially) chemi- 
cally differentiated structure. Third, there is an accretion phase 
during which the white dwarf grows in mass and increases its 
central density. At the accretion rates necessary to allow a SNIa 
explosion, the duration of this last phase can be much less than 
1 Gyr. For nearby SNIa, which presumably sample a range of 
progenitor ages, the sedimentation of 22 Ne during the second 
phase might introduce a larger variability of luminosity than for 
SNIa exploding at higher redshift, for instance z ~ 1. 

From the theoretical point of view, there is a long his- 
tory of calculations that have tried to unravel the influence of 
22 Ne s edimentation on the ph ysics of SNIa (e.g. iBravo et aD 
119921) . iBildsten & Halll d200lh proposed to use the production 
of 54 Fe in SNIa as a n indirect test of the sedimentation of 22 Ne. 
iPiro & Ch ang (2008) analyzed the impact of 22 Ne sedimentation 
on the size of the convective zone during the last stages of carbon 
simmerin g, which set the initi al conditions for a SNIa explosion 
(see also lPiro & Bild sten 2008). They claimed that to have an ap- 
preciable effect it would be necessary for all the 22 Ne to diffuse 
into the convective core prior to carbon ignition. In t heir mod- 
eling o f sub-Chandrasekhar white dwarf explosions, ISim et all 
(1201 Oh had to resort to a hypothetical gravitational settling of 
22 Ne to reproduce the central concentration of neutron-rich iso- 
top es inferred from ob serv ations of SNIa. 

iGarcia-Berro et all d2008l). lAlthaus et alj d2010h . 
iRenedo et al] d2010h . and Garcfa-Berro et al] i 20101) com- 



puted realistic white dwarf cooling sequences with updated 
physics, in which the evolution was followed simultaneously 
with the several possible chemical differentiation processes: 
gravitational settling of 22 Ne during the liquid phase and chemi- 
cal separation during the liquid-solid phase transition. It is, thus, 
interesting to know the consequences of these state-of-the-art 
white dwarf evolutionary calculations on the outcome of the 
thermonuclear explosions of these objects, to see whether they 
meet the expectations raised by the above-mentioned models of 
SNIa. The plan of the paper is as follows. In the next section, 
we briefly outline the chemical differentiation processes that 
take place during white dwarf cooling. The chemical profiles 
obtained at both the beginning and the end of the cooling pro- 
cess are used as input models for the thermonuclear supernova 
code. In Section [3] we report the results of the explosion of 
a Chandrasekhar-mass white dwarf, and the sensitivity of the 
supernova properties to the age of the progenitor. Section [4] is 
devoted to evaluating the impact of the chemical profiles on 
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Fig. 1. Chemical profile of the white dwarf at the beginning 
of the cooling process (red) and the the end of crystallization 
(green). The three abundances shown for each model belong 
(from top to bottom at M = 0) to: I6 0, I2 C, and 22 Ne. The mass 
fraction of 12 C in the outermost layer is 0.61 (red line) and 0.86 
(green line). 



sub-Chandrasekhar white dwarf explosions, where we account 
for their possible impact on the formation of a carbon deto- 
nation at the top of the carbon-oxygen core as a consequence 
of a helium detonation (the double-detonation scenario for 
sub-Chandrasekhar white dwarfs). Finally, in Section [5] we 
summarize our major findings and draw our conclusions. 

2. Chemical differentiation during white dwarf 
cooling 

The white dwarf evolutionary models used in this work were 
extracted from the full evolution of a 1.0 Mq cooling se- 
quence computed w ith the LPC0DE stellar evolutionary code (see 
Altha us et al.ll20 1 Ol for details). Particular attention was devoted 
to the treatment of the abundance changes resulting from the var- 
ious physical processes acting during the cooling phase, namely 
element diffusion in the outer layers, carbon-oxygen phase sep- 
aration upon core crystallization, and the slow 22 Ne sedimenta- 
tion in the liquid regions. The last process is particularly rele- 
vant to the present work. 22 Ne sedimentation towards the center 
of the white dwarf results from an imbalance between gravita- 
tional and electrical forces caused by their being tw o extra neu- 
trons in the 22 Ne nucleus relative to A, _ 2Z, nuclei dBrayo et al.l 
[l99llDel ove & Bildste nl2002tlGarcfa-Berro et al.l2008l) . In our 
simulations, time-dependent 22 Ne sedimentation was treated in 
a self-consistent w ay with the white d warf evolution, and we 
refer the reader to lAlthaus et al.l d2010l) for details. In partic- 
ular, for the liquid regions, we adopt the diffusion coefficient 
D = 7.3 x 10- 7 T/g l/2 Z F 1/3 cm 2 s 1 (see iDelove & Bildsten! 
2002). The chemical stratification of our starting white dwarf 
configuration at the beginning of the cooling track (depicted in 
Fig. [TJ is the result of the complete evolution of an initially 
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Fig. 2. Same as Fig.Q]but f° r a diffusion coefficient artificially 
increased by factor of five. 



5.0 M0 model star evolved from the zero-age main sequence, 
through the core hydrogen-burning phase, the helium burning 
phase, and the thermally pulsin g asymptotic gian t branch phase, 



Renedoetal 



2010). The initial 



to the white dwarf stage (see 

mass abundance of 22 Ne in the core is 0.02 (by mass). 

Noticeable changes in the chemical abundance distribution 
take place as the white dwarf evolves. During the evolutionary 
stages prior to core crystallization, abundance changes in the 
white dwar f core result from 2 2 Ne se dime ntation. As demon- 
strated by iDelove & Bildstenl (12002b and iGarcfa-Berro et al.l 
(2008), 22 Ne sedimentation is a very slow process that influences 
the evolution of white dwarfs only after long enough times have 
elapsed. During the evolutionary stages where most of the white 
dwarf remains in a liquid state, this process causes a strong de- 
pletion of 22 Ne in the outer region of the core, and an enhance- 
ment of its abundance in the central regions of the star. This can 
be noted from Fig. [TJ When most of the white dwarf has crystal- 
lized — at « 6 Gyr — there remains no trace of 22 Ne in the outer 
parts of the core. This is more noticeable in the case of a more 
efficient diffusion, as shown in Fig. [2] If we multiply the (rather 
uncertain) diffusion coefficient by a factor of 5, we indeed find 
that 22 Ne diffuses so deep into the core that no trace of this ele- 
ment is found in layers even as deep as 0.4 Mq below the stellar 
surface. 

Finally, we note that the carbon-oxygen distribution becomes 
strongly modified by the time the star has ended its crystalliza- 
tion process. This is a result of the carbon-oxygen phase separa- 
tion and the induced mixing episodes i n the outer liquid l ayers 
that take place upon crystallization (see lRenedo et al.ll2.Q10h . In 
particular, during the crystallization process, the oxygen abun- 
dance in the crystallizing region increases, and the overlying liq- 
uid mantle becomes carbon-enriched as a result of mixing in- 
duced by a Rayleigh-Taylor instability at the region above the 
crystallized core. Once the crystallization process is completed, 
the chemical profiles of all the elements become frozen, includ- 
ing that of 22 Ne. 



3. Delayed detonations of Chandrasekhar-mass 
white dwarfs 

Here we report the results of the explosion of a massive white 
dwarf using the chemical profiles obtained at the beginning and 
the end of the cooling process. The pre-supernova model is a 
cool, isothermal, Chandrasekhar-mass white dwarf built in hy- 
drostatic equilibrium with a central density p c = 3 x 10 9 g cm" 3 . 
The central 1 Mq has the same chemical composition as the 
carbon-oxygen core obtained in the cooling sequences, while the 
composition of the envelope of mass Mch - 1 Mq is the same 
as that of the outermost shell of the carbon-oxygen core. As we 
wish to maximize the impact of the chemical profiles on the su- 
pernova outcome, we do not take into account here convective 
mixing during the pre-supernova carbon simmering phase. 

The present m odels are based on the delay ed-detonation 
(DDT) paradigm dKhokhlovl Il99ll) . in which thermonuclear 
combustion initially proceeds through a subsonic deflagration 
until it makes a transition to a supersonic detonation wave. 
The location of the deflagration-detonation transition is usually 
parametrized by its density, pddt- The supernova hydrodynam- 
ics and nuc leosynthetic ( post-p roces sing) codes we use ar e the 
same as in iBravo et all d 19961) and iBadenes et al.l (|2003). We 
sketch here, for completeness, the method used by the hydrody- 
namics code to simulate DDT models and calculate the nuclear 
energy generation rate. In the DDT models, the flame propa- 
gates initially as a deflagration, with a velocity fixed at a con- 
stant fraction, 3%, of the local sound velocity. When the flame 
density ahead of the flame reaches the prescribed transition den- 
sity, pddt, the flame front is accelerated artificially to a large 
fraction of the sound speed, resulting in the subsequent forma- 
tion of a detonation. Changes in the chemical composition are 
followed during the propagation of the burning front in the de- 
flagration and detonation modes using an a-network from 4 He 
to 28 Si plus the conversion of 28 Si to 56 Ni in a single step. When 
the temperature of a mass shell exceeds 5.5 x 10 9 K, nuclear 
statistical equilibrium (NSE) is assumed. Once attained, NSE is 
maintained as long as the temperature remains above 2 x 10 9 K, 
providing a nuclear energy generation rate accurate enough for 
the hydrodynamical simulations. Weak interactions during NSE 
determine the evolution of the electron mole number, Y e , 



"dT 



(1) 



where Xu accounting for all kind of weak interactions, and 
the molar fractions, Y,, are set by the NSE equations. The 
final nucleosynthesis was computed separately with th e nu- 
clear reaction network described by Brav o et al.l 0993), with 
updat ed reaction rates taken f rom the REACLIB compilation 
( Raus cher & Thielem ann 2000), using the temperature and den- 
sities of each mass zone provided by the hydrocode. 

Bolometric light curves were obta ined by means of the code 
described in IBravo et all (1 1 9931 d?96). In general, the bolomet- 
ric light curves obtained with this code during the pre-maximum 
and maximum phases are in fairly good agreement with those 
comp uted by directly s olving the radiative transfer equations 
(see lHoflicheTaIlll993h . 

Table[TJgives details of the models we have computed so far. 
Each model is characterized by two parameters: its cooling time 
(age), and the deflagration-detonation transition density, pddt- 
This density is the main unknown in the DDT SNIa models. In 
one-dimensional calculations such as those reported here, pddt 
is usually a free parameter. While multi-dimensional SNIa mod- 
els may be able to emulate the physical dependences of pddt, in 
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Table 1. Results of Chandrasekhar-mass explosions. 



Age" 


Pddt 


M( 56 Ni) 


M(lME) b 


K c 


d 


(Gyr) 


(g cm j ) 


(Mq) 


(Mo) 


(10 51 erg) 


(mag) 


0.6 


1.6 x 10' 


0.40 


0.67 


1.26 


-18.78 


8.0 


1.6 x 10 7 


0.43 


0.70 


1.37 


-18.84 


1.0" 


1.6 x 10 7 


0.41 


0.67 


1.27 


-18.79 


9.6 e 


1.6 x 10 7 


0.43 


0.69 


1.38 


-18.88 


1.0" 


3.0 x 10 7 


0.71 


0.42 


1.42 


-19.45 


9.6 e 


3.0 x 10 7 


0.74 


0.42 


1.51 


-19.50 


0.6 


2.4 x 10 7 / 


0.59 


0.52 


1.37 


-19.23 


8.0 


3.7 x 10 7 « 


0.83 


0.34 


1.54 


-19.63 



Notes. All models: central density, p t . = 3.0 x 10 9 g cirr 3 . 

<a) Time since the beginning of the cooling sequence. 

{b) Synthesized mass of intermediate-mass elements, from Si to Ca. 

{c} Kinetic energy of the ejecta. 

{d) Bolometric magnitude at maximum. 

(e) Model with the diffusion coefficient of 22 Ne increased by a factor 5. 
W) Deflagration-to-detonation transition density estimated from the 
chemical composition at the flame front (Eq. [2^: X( I2 C) = 0.36, 
X( 22 Ne) = 0.020. 

te) Deflagration-to-detonation transition density estimated from the 
chemical composition at the flame front using Eq. {2j: X( 12 C) = 0.24, 
X( 22 Ne) = 0.022. 



principle, state-of-the-art models are still nowada ys not able to 
do it ( for a recent attempt in this direction see, e.g. lJackson et alJ 
120101) . Thus, we repeated the simulations for several values of 
Pddt, as shown in Table [TJ In some calculations, we used ini- 
tial models resulting from using a diffusion coefficient of 22 Ne 
during white dwarf cooling artificially increased by a factor of 
five. For each combination of pddt and diffusion coefficient, we 
computed two models, one belonging to the beginning of crys- 
tallization of the white dwarf ("young" model) and another at an 
advanced state of crystallization (~ 95 - 96%, "old" model). 



3.1. Detonating at a fixed p DD j 

The first two rows of Table [TJ give the mass of radioactive 56 Ni, 
the mass of intermediate-mass elements (IME, mostly composed 
of Si and S), the kinetic energy of the ejecta, and the bolomet- 
ric magnitude of the supernova at maximum, for two models 
computed at the beginning and the end of the white dwarf cool- 
ing sequence, with pddt = 1-6 x 10 7 g cm 4 . The results are 
quite insensitive to the age of the pre-supernova, the difference 
in brightnesses being ~ 0.06 mag. The differences in both the 
yield of IMEs and the kinetic energy are correspondingly small, 
on the order of 5% and 9%, respectively. 

3.1 .1 . Enhanced diffusion of 22 Ne 

The third and fourth rows of Table Q] give the results of models 
identical to the ones just discussed but for the enhanced diffu- 
sion coefficient of 22 Ne during white dwarf cooling. The results 
are nearly identical to those obtained with the nominal diffu- 
sion coefficient. Thus, any uncertainty in the diffusion coeffi- 
cient does not translate into appreciable differences in the ob- 
servational properties of SNIa. These conclusions are indepen- 
dent of the value of pddt used in the supernova calculations, as 
can be seen in the next two rows where we show the results for 
Pddt = 3.0 x 10 7 g cirr 3 . 

Figures [3] and [4] show the final chemical composition of the 
ejecta as a function of the ejecta mass (as a Lagrangian coordi- 




0.5 

Lagrangian mass (M ) 



Fig. 3. Final chemical composition of the ejecta for the model 
detonated at the beginning of the cooling process with a fixed 
Pddt = 1.6 x 10 7 g cirr 3 and enhanced diffusion. The abun- 
dances shown reflect the elemental composition after radioactive 
disintegrations. See the online edition of the journal for a color 
version of the figure. 



nate), for the two models with enhanced diffusion coefficient and 
Pddt = 1 .6 x 10 7 g cm 4 . The main difference between the initial 
models (see Fig. [2j lies in the concentration of 22 Ne in the cen- 
tral 0.6 M0, which leads to a larger 22 Ne mass fraction (~ 50%) 
in the oldest model. However, the chemical profiles of the ejecta 
of both models are nearly the same below this Lagrangian mass. 
Between ~ 0.6 Mq and 1 M©, the oldest model (8.0 Gyr age) 
has slightly larger mass fractions of Si and S and less Fe and 
Cr, although the differences are quite modest. The reason for 
these differences is the following. As the thermonuclear wave 
propagates through the central regions of the white dwarf, it 
completely incinerates the material up to a composition in nu- 
clear statistical equilibrium. The oldest model having a larger 
central neutron excess (larger 22 Ne mass fraction) releases less 
energy during the incineration of carbon and oxygen in reaching 
NSE, which causes a slightly slower expansion during the ini- 
tial phases of the explosion, and higher densities at the position 
of the flame. As a consequence, the subsequent combustion of 
the layers above a ~ 0.6 Mq Lagrangian mass is able to reach a 
more advanced combustion stage than if the expansion had been 
faster (as in the youngest model, of 0.6 Gyr age), and more Fe 
and less IMEs are produced in these layers. 

3.2. Transition density function of the chemical composition 

We now speculate wh ether pddt is a funct i on of the local 
chemical composition Jchamulak et al.l 120071; IWooslevI [2007; 
iBravo et"al ] l2010tlJackson et al.l l2010fl In considering this sce- 



1 The models discussed in this section are based on the nominal value 
of the diffusion coefficient. 
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0.5 

Lagrangian mass (M Q ) 



0.5 

Lagrangian mass (M ) 



Fig. 4. Same as Fig. [3] but at the end of the crystallization phase. 



nario, we scaled the transitio n density as a func tion of the local 
chemical composition as (see lBravo et al.ll2~010i for details) 



Pddt oc X( 12 Cy h3 (1 + 129/7)- - 6 



(2) 



where the neutron excess, rj — 1 —2Y e , is related to the 22 Ne mass 
fraction by 



tj = 



X{ zl Ne) 

il 



(3) 



Equation © implies that a reduction in the carbon mass frac- 
tion translates into a higher pddt, while an increase in the 22 Ne 
mass fraction has the opposite effect. The results we obtained 
are given in the last two rows of Table [1] and in Figs. [5] and |6l 
where the final chemical composition of the ejecta is shown for 
the models corresponding to both the beginning and the end of 
the white dwarf cooling process. 

Allowing pddt to be a function of the local (at the position of 
the flame front) chemical composition makes a great difference. 
The transition from deflagration to detonation takes place at a 
Lagrangian mass ~ 0.22 IVLtJI The abundances of 12 C and 22 Ne 
at this location are 0.36 and 0.020, respectively, for the youngest 
model, and 0.24 and 0.022 for the oldest model (see Fig. Q]). 
The main effect comes from the reduction in the carbon mass 
fraction with age, which causes an increase of ~ 63% in pddt- 
The modest variation in the neon abundance at the flame front 
only makes the transition density decrease by ~ 1%. Finally, the 
difference in pddt as a function of the progenitor age translates 
into a SNIa brightness difference of 0.40 mag. 

In these calculations, the production of IMEs by SNIa is 
strongly affected by the age of the progenitor, since it is reduced 
by a factor of ~ 0.65 for the oldest model. As can be seen in 



2 The location of the DDT can be identified in Figs. [3] to [6] with the 
sudden drop in the Ni profile slightly above 0.2 Mq. This discontinuity 
is a numerical artifact of the procedure used to transmute a deflagration 
front into a detonation wave. 



Fig. 5. Chemical composition of the ejecta for the model with 
Pddt function of the mass fractions of 12 C and 22 Ne as given by 
Eq. (f2]i at the beginning of the cooling process: X( l2 C) - 0.36, 
X( 22 Ne) = 0.020 (see Table Q). 




n mass (M ) 



Fig. 6. Chemical composition of the ejecta for the model with 
Pddt function of the mass fractions of 12 C and 22 Ne as given by 
Eq. (f2]i at the end of the crystallization phase: X( l2 C) = 0.24, 
X( 22 Ne) = 0.022 (see Table [T). 



Figs. |5] and [6] the larger pddt is, the more extended the Fe pro- 
duction region and the more squeezed the region rich in IMEs 
is. 
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4. Explosion of sub-Chandrasekhar white dwarfs 

The sub-Chandrasekhar model of SNIa 



Table 2. Results of sub-Chandrasekhar explosions. 



1980 


; iNomoto 


1982; Wooslev et alj 


1986; Iben&Tutukov 


1991 


; Limonai & Tornambe 


1991; Livne & Glasnerl 


1991; 
chal- 


Wooslev & Weaver 1994; Garcfa-Senz et al. 1999) was 



lenged 14 years ago because its spectra was excessively blu e 
compared to SNIa observations dHoflich & Khokhlovl Q996). 
This disagreement had its origin in the assumption that a thick 
helium layer was necessary to detonate a sub-Chandrasekhar 
white dwarf. However, there has been renewed interest in this 
kind of explosion because of the possibility that they produce a 
new class of underluminous supernovae (Bild stenetalJ 120071) . 
and the finding that a white dwarf detonation is p ossible for 
lower helium shell masses than previously thought dFink et alJ 
120101) . This last p oint allows a better agreement with SNIa 
optical properties dKromer et al .] 1201 01) but the models still 
show some deficiencies, particularly in terms of their lack of a 
central region rich in stable Fe-grou p elements, as demanded 



by infrared observations of SNIa dHoflich et alJ I2004j). One 
possib le solution for this deficiency, proposed by Kramer et al. 
d201Q ). might be a higher neutron excess in the central regions 
of the white dwarf caused by gravitational settling of 22 Ne 
during white dwarf cooling. We test this hypothesis in the 
following. 

The initial models used in our sub-Chandrasekhar explo- 
sion calculations are composed of a carbon-oxygen core and 
a helium envelope. The core has the same chemical structure 
as the carbon-oxygen core in the 1 Mq white dwarf evolved 
with the LPCODE code (Sect. [2]), taken at different ages since the 
beginning of the cooling process. The envelope has a mass of 
0.080 Mq and is composed of 100% He. The mechanical struc- 
ture of the white dwarf that resulted from the cooling sequences 
was readjusted to reestablish the hydrostatic equilibrium tak- 
ing into account the envelope. As a result, the central density 
changed from ~ 3.5 x 10 7 g crrr 3 , as given by the cooling se- 
quences, to ~ 5.8 x 10 7 g crrr 3 in our initial models. The density 
at the base of the He envelope was ~ 2.1 x 10 6 g cm 4 . 

The models we computed are summarized in Table [2] The 
models were calculated at two different ages of the progenitor 
white dwarf, and the initial He detonation was started at two 
different locations: just at the base of the helium envelope, and 
at an altitude of 70 km above this. We also computed a reference 
model in which the composition of the core is uniform and has 
equal mass fractions of carbon and oxygen (first row in Table|2]). 

The precise point at which He detonates is not known with 
precision because of unresolved convective motions during the 
first stages of He burning at the base of the envelope, so the 
above procedure allows us to evaluate the impact of the chemi- 
cal differentiation of the white dwarf in two extreme conditions. 
Detonating helium at different altitudes can have consequences 
for the outcome of th e inward shock wave t hat is launched into 
the white dwarf core dLivne & Glasnerifl99ll) . If He ignition oc- 
curs at a high enough altitude above the core, the He detonation 
wave has time to develop well before it arrives at the core and it 
can be strong enough then to directly induce an inwardly mov- 
ing carbon detonation. This possibility depends as well on the 
carbon mass fraction at the outermost part of the core, which in 
turn depends on the age of the progenitor white dwarf. 



4.1. Edge-lit detonation of carbon 

We explored the conditions for the formation of a stable self- 
sustained carbon detonation induced by a He detonation, as 



Age" 




M( 5ft Ni) 


M(Fe, sta) c 


M(IME) 


K 


(Gyr) 


(km) 


(Mo) 


(Mo) 


(Mo) 


(10 51 erg) 


j 





0.70 


0.055 


0.22 


1.39 


0.6 





0.69 


0.054 


0.23 


1.34 


8.0 





0.68 


0.058 


0.24 


1.35 


0.6 


70 


0.63 


0.054 


0.28 


1.31 


8.0 


70 


0.64 


0.057 


0.30 


1.34 



Notes. All models: Mass of the carbon-oxygen core, M core = 1 Mq; 
mass of the He envelope M cnv = 0.080 Mq; central density, 
p c = 5.8 x 10 7 g cirr 3 ; density at the base of the He envelope, 
ptase = 2.1 x 10 6 gem" 3 . 

(a> Time since the beginning of the cooling sequence. 
{ - b) Altitude at which the He detonation starts. 
{c) Ejected mass of stable Fe-group nuclei. 

(<fl Reference model, with a homogeneous composition of the core of 
49% carbon and oxygen plus 2% 22 Ne. 



a function of the carbon mass fraction and the altitude at 
which He is ignited. A similar analysis was performed by 
iGarcfa-Senz et al.| (fl999). but they explored a range of 12 C mass 
fractions below 0.5, while those we are interested in are much 
larger (X( l2 C) - 0.61 in our young prog enitor model, and 
X( l2 C ) — 0.86 in the old progenitor model). IGarcfa- Senz et al. 
d 19991) found that the probability of a self-sustained detonation 
occurring in the core is very sensitive to the carbon mass frac- 
tion. 

In the present numer ical experiments, we applied the same 
methodology adopted in iBravo & Garcfa-Senzl d2009l and ref- 
erences therein). We followed the hydrodynamical and nuclear 
evolution of an isothermal uniform density sphere consisting of a 
central ball made of 99% 4 He and 1 % I2 C, surrounded by a much 
higher mass of the same chemical composition as the outermost 
layer of the carbon-oxygen core of the white dwarf, produced by 
the cooling sequences. The radius of the central ball plays the 
role of the altitude at which He ignites in the sub-Chandrasekhar 
models. The He detonation was started by incinerating a central 
region containing between 1% and 8% of the mass of the central 
helium-rich ball. 

The results of these numerical experiments are summarized 
in Table f3] where we have investigated the induction of carbon 
detonation at two densities: 2 x 10 6 g crrr 3 , which matches the 
density at the base of the helium envelope in our 1.080 Mq 
sub-Chandrasekhar models, and 6 x 10 6 g crrr 3 , which would 
be representative of models with a higher total mass. Here, 
a "weak detonation" means a carbon detonation in which the 
temperature is no t high enough to burn oxygen directly (see 
IGarcfa-Senz et al.1 1 19991) . In a sub-Chandrasekhar explosion, 
such an edge-lit weak detonation will strengthen as it moves into 
denser regions, reaching increasingly higher temperatures and 
allowing the completion of burning sequences from explosive 
oxygen burning until NSE. The successful launching of a carbon 
detonation wave is found to depend strongly on the density and 
the altitude of He ignition, and weakly on the carbon mass frac- 
tion. For instance, for a density of 2 x 10 6 g crrr 3 , the minimum 
altitude at which a carbon detonation develops in our simulations 
is 200 km for a composition of 50%:50% carbon and oxygen, 
and only 100 km for the two compositions richer in carbon that 
we explored. For the highest density we used, 6 x 10 6 g cm -3 , 
and an altitude of 5 km, a carbon detonation was obtained only 
for the largest carbon mass fraction, X( l2 C) = 0.81. Thus, the 
age of the progenitor system can determine the kind of dynam- 
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ical event following He ignition in a sub-Chandrasekhar white 
dwarf, i.e. either a central carbon detonation or a shell carbon 
detonation. 

4.2. Explosion models 

The effects of white dwarf crystallization and 22 Ne sedimenta- 
tion on sub-Chandrasekhar explosions can be seen in Table [2] 
First, comparing the three models that ignited He at the base of 
the envelope (first three rows, with h{ g = 0) we see that the im- 
pact of chemical differentiation is minimum, on the order of a 
few percent in 56 Ni and IMEs yields and in the kinetic energy. 
Second, comparing the two models with h[ g = 70 km, we ob- 
serve the same lack of dependence of the outcome of the explo- 
sion on the age of the progenitor system. 

A further comparison can be established between a young 
model in which the edge-lit detonation of carbon fails (and, later, 
carbon detonates at the center of the white dwarf), and an old 
model in which the edge-lit detonation of carbon is success- 
ful. Figures [7] (young model, second row in Table |2]l and [8] (old 
model, last row in Tabled show the evolution of the white dwarf 
in both cases. Even though the history of the explosion is quite 
sensitive to the location of carbon detonation, the final result is 
not so different. In the last panel of Figs. [7] and [8] we show the 
final mechanical, thermal, and velocity profiles of the ejecta, to- 
gether with the profile of the mass fraction of stable Fe-group 
elements (red dotted line). The concentration of stable Fe-group 
elements is clearly insensitive to the age of the progenitor system 
of the exploding sub-Chandrasekhar white dwarf. 

5. Conclusions 

We have investigated the sensitivity of SNIa energetics and nu- 
cleosynthesis to realistic chemic al profiles based on w hite dwarf 
cooling sequences calculated by Altha uset al.l (12.010b . The cool- 
ing sequences provide a link between the white dwarf chemi- 
cal structure and the age of the supernova progenitor system. 
Neither gravitational settling of 22 Ne nor chemical differentia- 
Table 3. Induced detonation of carbon. 



p 


X( l2 Cf 


Altitude" 


Detonation? 


(g cm" 3 ) 




(km) 




6x 10 b 


0.50 


5 


no 


6x 10" 


0.61 


5 


no 


6x 10 6 


0.86 


5 


yes 


6x 10" 


0.50 


10 


weak 


6x 10 6 


0.61 


10 


yes 


6x 10" 


0.86 


10 


yes 


6x 10 6 


0.50 


20 


yes 


2x 10" 


0.50 


50 


no 


2x 10" 


0.61 


50 


no 


2x 10" 


0.86 


50 


no 


2x 10 6 


0.50 


100 


no 


2x 10" 


0.61 


100 


weak 


2x 10 6 


0.86 


100 


weak 


2x 10" 


0.50 


200 


weak 


2x 10 6 


0.61 


200 


weak 


2x 10 6 


0.86 


200 


weak 



Notes. <n> Carbon mass fraction surrounding the central He ball. In all 
calculations, the oxygen mass fraction in this region was set to X( l6 0) = 
l-X(' 2 C). 

{b) Radius of the central ball made of 99% 4 He and 1% 12 C. 
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Fig. 7. Development of the detonation of a sub-Chandrasekhar 
mass white dwarf that explodes at the beginning of the cool- 
ing process. From left to right and top to bottom, the first three 
panels show the profiles of temperature in units of 10 9 K (short- 
dashed line), density in units of 10 7 g cm" 3 (long-dashed line), 
and velocity in units of 10,000 km s" 1 (solid line), at three times 
since the beginning of He detonation. In this model, the He det- 
onated at the base of the envelope, so it was unable to produce 
an edge-lit detonation of carbon (see Tables [2] and [3]). Instead, 
an initially weak inwardly moving shock traveled through the 
white dwarf (second panel). Shortly before reaching the center, 
the shock strengthened enough to produce a carbon detonation 
that traveled outwards processing most of the matter on top (third 
panel). The last panel shows the final profiles of log(T) (short- 
dashed line), log(10 5 p) (long-dashed line), velocity (solid line), 
and 50X(Fe, stable), where X(Fe, stable) is the mass fraction of 
stable Fe-group nuclei (red dotted line). The discontinuity in the 
thermal profile in the last panel is due to the local deposition of 
the photons emitted in the radioactive disintegration of 56 Ni. 



tion of I2 C and I6 during white dwarf crystallization have a 
sizable impact on the properties of SNIa, unless there is a direct 
dependence of the flame properties on chemical abundances. If 
the density of transition from deflagration to detonation in SNIa 
did not depend on the chemical composition, the variation in the 
supernova magnitude with age produced by chemical differenti- 
ation would be as small as ~ 0.06 magnitudes for an age differ- 
ence of ~ 7.4 Gyr. We emphasize that these results have been 
obtained by neglecting mixing during the pre-supernova carbon 
simmering phase. If this mixing process is efficient, it will erase 
any trace of the chemical separation achieved during white dwarf 
cooling, thus the chemical separation will leave no imprint on the 
supernova properties. Our results therefore represent the maxi- 
mum possible effect that can be expected from the gravitational 
settling of 22 Ne and ls O. 

If the density of transition from deflagration to detonation 
in SNIa is a function of the chemical composition, and if we 
neglect convective mixing during the pre-supernova carbon sim- 
mering phase, the difference in maximum magnitude between a 
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Fig. 8. Development of the detonation of a sub-Chandrasekhar 
mass white dwarf that explodes at the end of the crystalliza- 
tion phase. The meaning of the different curves is the same as 
in Fig. [7] In this model, the He detonated at an altitude of 70 km 
from the base of the envelope, and launched a second inwardly 
moving carbon detonation (see Tables [2] and [3]), as can be seen 
already in the second panel. The carbon detonation processed all 
the matter down to the center of the white dwarf and, thereafter, 
sent a weak shock wave outwards (visible in the third panel at a 
Lagrangian mass of ~ 0.5 Mq). 



white dwarf exploding at the beginning of the cooling process 
and another one that has experienced substantial crystallization 
can be as large as ~ 0.4 mag. For this variation in magnitude 
to take place, it suffices that the central part of the white dwarf 
has crystallized (22% in mass in our models), so that the transi- 
tion from deflagration to detonation takes place in a region that 
has been depleted in carbon. For a 1 Mq white dwarf, the crys- 
tallization of the central 22% in mass occurs at an age of only 
1.25 Gyr. 

The physics of 22 Ne sedimentation is nowadays sufficiently 
well-known that the calculations presented in this paper cannot 
experience much variation. We have also explored the impact of 
using an unrealistically large diffusion coefficient, which nev- 
ertheless had a very small impact on the observable properties 
of SNIa explosions. We stress that to obtain a large concentra- 
tion of 22 Ne in the central layers of a white dwarf the star would 
have to be made of pure carbon, which is not applicable in the 
more general case of carbon-oxygen white dwarfs. Thus, our re- 
sults concerning the lack of imprint of the 22 Ne sedimentation 
on SNIa are quite robust. 

In the light of the present results, the only possible ways in 
which chemical differentiation might affect the thermonuclear 
explosion of a white dwarf are either through a composition de- 
pendent density of transition from a deflagration to a detonation, 
or through composition dependent thermonuclear runaway con- 
ditions (usually taken to be the initial conditions for SNIa sim- 
ulations). Linking explosion models to the last stages of white 



dwarf evolution remains one of the main challenges of SNIa the- 
ory. 

Finally, we highlight the importance of developing realistic 
pre-supernova evolutionary calculations to understand the con- 
straints imposed by observations on SNIa models. In the case 
of sub-Chandrasekhar explosions, we have proven that the sed- 
imentation of 22 Ne does not efficiently increase the neutron ex- 
cess in the central region of the white dwarf prior to the explo- 
sive event. Thus, this sedimentation cannot be responsible for 
the production of the stable Fe-rich core inferred from observa- 
tions of SNIa. However, the crystallization of the white dwarf 
entails an increase in the carbon abundance in the outermost lay- 
ers of the carbon-oxygen white dwarf core, which may have con- 
sequences for the formation of a double detonation (outwards in 
He-rich matter, inwards in the C-O core) following helium igni- 
tion close to the base of the helium mantle. 
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